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Abstract. Compact object mergers are one of the fa- 
vorite models of GRBs. It has been noted that in opposi- 
tion to the coUapsars, compact object mergers do not nec- 
essarily take place in the host galaxies, and may travel out- 
side of them. With the discovery of afterglows and identi- 
fication of host galaxies one can measure the distribution 
of GRBs with respect to their host galaxies. This distri- 
bution has been calculated using different population syn 
thesis codes , and for different ga l actic potentials (Bloom 
et al, 199g| ; [Bulik et al., 1999a| ; [Fryer et al., 1999| ). In 



this paper we compare the distributions of different types 
of compact object binaries: double neutron star systems 
(NS-NS), black hole -neutron star systems (BH-NS) and 
double black holes (BH-BH). We calculate the orbits and 
distributions of the projected distances on the sky for two 
extreme cases: a massive galaxy like the Milky Way, and 
empty space (corresponding to e.g. a globular cluster), 
and consider a wide range of possible kick velocity distri- 
butions. We find that BH-NS are more likely gamma-ray 
burst counterparts, since they lie close to the host galax- 
ies, contrary to the NS-NS binaries. 
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1. Introduction 



gamma rays: 



The discovery of X-ray (Costa et al., 1997) and later opti- 
cal afterglows (Groot et al., 19971) of gamma-ray bursts 
led to jdeiitification of p-armna-ray burst host fyalaxies 



associations are not obvious. Compact object binaries may 
live for quite a long time before they merge, and given the 
possibility that they could have high velocities they may 
travel away from the place they were formed. 

The distribution of compact object mergers can be 
found using stellar population codes. The main problem 
of this approach is that such codes contain a number of 
poorly known parameters, which may affect the results. 
One of the most important parameter is the kick velocity 
a newly born compact object receives at birth. 

The distribution of double neutron star systems 
around galaxies has been calculated for a few types of 
galaxies by Bloom et al. (1999). Bulik et al. (1999a) calcu- 
lated such distribution for the case of a massive galaxy like 
the Milky Way, and for the case of empty space, using four 
different kick velocity distributions. These studies consid- 
ered only binaries containing neutron stars. We used an 



assumption (Bulik ct al., 1999a; Bulik et al., 1999b) that 



all supernovae lead to formation of a 1.4 Mq neutron star. 
However, mergers of binaries containing a black hole are 
now more favored for GRBs. One reason for this is en- 
ergetics, GRB990123 had an equivalent isotropic energy 
release of 10*^^ ergs. The energetic requirements go down 
when considering that the relativistic outflow from the 
central fireball is not isotropic but beamed. On the other 
hand not all kinetic energy in the outflow can be converted 
into gamma rays, and therefore the energetic requirement 
for the GRB central engine will go up. Therefore it is im- 
portant to investigate the distribution around galaxies of 
binaries containing black holes. 



In this paper we extend the results of (Bulik et al 



( prootlet a1 1 997a[) Such host galaxies were expected m 1999a[ ), to include the case of compact object binaries con- 



the cosmological model, since most of the theories identi 
fied GRBs with some, perhaps extreme, stages of stellar 
evolution. However, various physical models of GRBs gave 
different predictions regarding the location of GRBs with 
respect to the host galaxies. In the framework of the col- 
lapsar model, and in general all models that relate GRBs 
to the final stages of evolution of massive stars ( polgate" 
196S), one expects that GRBs are found in the star form- 
ing regions. In the second class of models where bursts are 
associated with mergers of compact object binaries such 



taining black holes. In section 2 we describe the model for 
population synthesis and galactic potential used in this 
paper, in section 3 we present the results, and we summa- 
rize this work in section 4. 

2. The model 

2.1. Population synthesis code 

We use the populatio n syn thesis code described in detail in 
Belczyiiski & Bulik (1999). Within this model we assume 
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that the distribution of the masses of the primary stars is 

with 10 Mq < M < 50 Mq. Such a range was chosen be- 
cause we only want to deal with stars that may undergo 
a supernova explosion. The lower limit {Msn — 10 Mq) 
might seem high if one believes that large convective core 
overshooting is required in order to reproduce observations 
of single stars (e.g. observed wide width of main sequence 
bands) . Large convective core overshooting stellar models 
will result in M5jv= 6-8 Mq. But following the arguments 
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of (Vanbeveren et al., 1998) we note that the effects like 
stellar rotation or enhanced stellar winds may also explain 
the observations of single stars without need of large con- 
vective core overshooting. If one uses small convective core 
overshooting then the limiting mass of ZAMS star to be- 
come supernova becomes larger than 9 Mq. Change of 
MsAT to a lower value would increase the number of NS- 
NS binaries as compared to the number of BH binaries. 
This would not change our results as in this paper we work 
only on distribution of a given class of mergers around its 
host galaxy. 

The distribution of the mass ratio q (secondary to pri- 
mary mass) is 

<i>(g) oc const 

with < <7 < 1. The distribution of the initial binary 
eccentricity e is S(e) — 2e with < e < 1, and the distri- 
bution of the initial semi-major axis a used in population 
synthesis codes is flat in the logarithm, i.e. T{a) oc 
with the maximum Umax — lO^i?0. 

Very little is known observationally about the kicks 
newly born black holes receive. While studies of pulsars 
indicate that substantial kicks are possible for newly born 
neutron stars, we believe that black hole kicks are smaller. 
The physical reason could be just that black holes are 
more massive, and also that asymmetric neutrino emission 
may not take place when a black hole is formed. Lipunov 
et al. ( |1997| ) use a Gaussian parameterization of the kick 
velocity and assume that for black holes it is proportional 
to the mass lost in the explosion. Fryer et al. ( 1999| ) show 
results for several cases, using the same functional form of 
the kick distribution for black holes and neutron star and 
assuming that black hole kicks are ten times smaller than 
those for neutron stars. 

We parameterize the distribution of the kick velocity 
a newly born compact object receives in a supernova ex- 
plosion by a three dimensional Gaussian with the width 
a„, and consider several values of this parameter. Such a where Mch — 1-4 Mq is the Chandrasekhar mass, and 



Fig. 1. The population of the compact object binaries 
in the space spanned by their masses. The three regions 
correspond to different initial mass ratios: the binaries 
with Ml reaching up to 10 Mq and M2 < Mi, origi- 
nate from systems with the initial mass ratio q < 0.88, 
and the binaries along the line Mi « M2 originally had 
nearly equal masses q > 0.95. The sparse population of 
systems for which M2 > Mi originates from binaries with 
0.88 < g < 0.95. 



2.2. Compact object formation in a Sn explosion 

In general, supernovae from massive stars should lead to 
formation of more massive objects, perhaps black holes, 
and supernovae from less massive stars to lighter com- 
pact objects like neutron stars. The boundary between 
these two possibilities and the mass of a compact object 



formed in an explosion is not clear. Lipunov et al. (1997) 
assumed that a neutron star with a mass 1.4 Mq is formed 
for progenitors with core masses below 35 Mq at the time 
of explosion, and for larger masses black holes are formed 
with the mass k^h x 35 Mq where fcf,/,. = 0.3. The pa- 
rameter fcf,/j is varied between 0.1 and 1 (Lipunov, 1997). 



Tout et al. ( |1997| ) used a different prescription: neutron 
stars are born with masses equal to Mch + -^pro/50. 



choice allows to compare the distributions of binaries con- 
taining black holes vs. the double neutron stars that could 
be formed with different kick velocities. 

For a detailed description of the population synthesis 
code we refer the reader to ( Belczynski and Bulik, 19"99| ). 
There is only one modification to the code, i.e. in the 
choice of the mass of a compact object formed in a super- 
nova explosion. This is described in the next subsection. 



Mpro is the progenitor mass at the time of collapse, and 
assume that object above 1.8 Mq are black holes. Fryer 



et al. (1999) presented numerical simulations of core col- 
lapse in massive stars, and found that neutron stars are 
formed for progenitors with initial masses below 22 Mq, 
and black holes for m asses above this va lue. In a popu- 
lation synthesis code ( Fryer et al., 199£ ) assume simply 
that a compact object with a mass \Mp, where Mp is the 
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progenitor mass at the time of collapse. On the other hand 
Ergma and van den Heuvel ( 1998| ) analyze the population 
of known X-ray binaries and argue that the initial mass of 
the progenitor is not the sole parameter determining the 
mass of a compact object formed in a supernova explosion, 
but probably other factors like rotation or the magnetic 
field may play an important role. 

Given such a range of uncertainties we use the follow- 
ing prescription for the compact object mass (either neu- 
tron star or black hole): the compact object mass is equal 
to half of the mass of the final helium core of progeni- 
tor star. This approach gives higher masses of compact 
objects as compared to the prescriptions discussed above. 
However, simulations of cor e collapse during supernova ex- 
plosions ( Fryer et al., 1999 ) show that the compact object 
might accrete back some of the material firstly ejected in 
explosion for a fraction of the range of stellar masses we 
are concerned in our work. Although the predictions how 

much mass is arcrptpc] in a given r3sp are still nnrertain 



and caVi not he easily parampteriypH hnt wp bplievp that 

some enhancement of compact object mass is reasonable. 

Standard equations of state for neutron stars give the 
maximum mass of neutron star to be between 1.9 and 
2.6 Mq, although the rotation may increase it by about 
another 0.2 Mq. We have used value of 2.4 MgUn as the 



maximum neutron star mass in our previous work (Bulik 



et al., 1999a). However, there are equations of state for 
neutron stars which go as high as 3 Mq or even up to 
3.2 Mq if rotation is included ( |Cook et al., 1994] ). Thus in 
this work we assume that the maximum neutron star mass 
is 3.0 Mq. This number was not crucial in the previous 
paper, as we treated both groups, namely double neutron 
stars and black hole neutron stars binaries together and in 



(Bulik et al., 1999a) we assumed that all supernovae lead 
to formation of a 1.4 Mq neutron star. In present work 
this number plays an important role and as we are mostly 
interested in black hole binaries we want to be sure that 
we do not include any neutron stars as black holes in the 
present calculations. Thus we use 3.0 Mq as the upper 
limit on neutron star mass and consequently the dividing 
line between neutron star and black hole. 

2. 3. Gravitational potentials 

We follow the approach we have used in our previous pa- 
per ( Bulik et al., 1999a ), and consider two extreme cases: 
propagation in a potential of a massive galaxy like the 
Milky Way, and propagation in empty space (correspond- 
ing to e.g. a globular cluster origin). The potential of a 
massive galaxy consists of three components: bulge, disk, 
and halo. To model the bulge and disk potentials we use 
the potential model proposed by (Miyamoto and Nagai 
1975|): 



$(i?,z) = 



GM, 



+ (a, + v/i^Tfe|)2 



(1) 



where and hi are the parameters, M is the mass, and 
R = ^ + y^- The dark matter halo potential is spheri- 
cally symmetric 



$(r) 



GMh 



- In 1 + — + — atan — 



where Tc is the core radius. The halo potential corresponds 
to a mass distribution p = Pc/[^ + (''/^c)^], and we intro- 
duce a cutoff radius rcut ~ 100 kpc beyond which the halo 
density falls to zero, in order to make the halo mass fi- 
nite and the halo gravitational potential is $(r) oc 
when r > Vcut- We use the following values of the pa- 
rameters derived for the Milky Way (we assume that the 
Milky Way is a good example of a massive galaxy; the 
bulge potential {i = 1): oi = Okpc, bi — 0.277 kpc, A/i = 
1.12 X IO^^A^q; the disk potential (i = 2): 02 = 4.2 kpc, 
62 = 0.198 kpc. Ah = 8.78 x 10^" Mq; the halo poten- 
tial: rr, = 6.0 kpc, and Mh ^ 5.0 x 10^° Mq ([Paczynski 



1990; Blaes and Rajagopal, 1991 ). The distribution of stel- 



lar initial positions in the model galaxy is a double ex- 
ponential P{R,z) — Rexp{—R/Rexp)exp{—z/Z(,xp), with 
Rexp = 4.5 kpc, Zexp = 75 pc, and we cut the distribu- 
tion at Rmax = 20 kpc dBulik et al., 1998| ). The initial 
velocity of a binary is assumed to be equal to the local ro- 
tation velocity in the galactic disk. After each supernova 
explosion we add the velocity the system received, which 
is calculated in the population synthesis code. We follow 
the trajectories of the binaries until they merge. 



3. Results 



In order to see the general properties of the population 
of compact object binaries we consider the case with no 
kick velocities ay = Okms~^. In Figure |l| we present the 
density of binaries as a function of their masses. The popu- 
lation in Figure]^ has three components. The main compo- 
nent are the binaries which originate from systems with 
small mass ratios q < 0.88, the binaries along the line 
Ml w M2 originally had nearly equal masses q > 0.95 
and the population of systems for which M2 > Mi origi- 
nates from binaries with intermediate mass ratios 0.88 < 
q < 0.95. These components come from different regimes 
of initial binary mass ratio, which to some degree sets 
the subsequent binary evolution, as discussed earlier by 
(Bethe and Brown, 1998). The detailed description of evo- 
lution in each of these regimes can be found in Belczynski 
& Buhk ( p^99| ). 

When considering higher kick velocities the basic shape 
of the diagram of Figure |l] does not change substan- 
tially. We draw lines corresponding to Mi — 3Mq and 
M2 = 3Mq, the upper limit on the maximal mass of a 
neutron star. We predict a large population of compact 
object binaries that contain a black hole. While we do not 
expect to see double black hole binaries, it is puzzling that 
no black hole neutron star binaries are known among bi- 
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Fig. 2. The distributions of binaries in the plane spanned by the center of mass velocity v and merger time tmrg- In 
each panel we also plot the following lines: the horizontal dashed line corresponding to the Hubble time (ISMyrs); two 
solid lines: the vertical corresponding to f = 200 km s"""^ - a lower limit on the typical escape velocity from a galaxy, 
and the line corresponding to a constant value of v x tmerge = 30kpc. The panels correspond to the kick velocity 
distributions width of 0km s~^, 200km s~^, and 800km s~^ from top to bottom respectively. In each panel we present 
the distribution of NS-NS, BH-NS and BH-BH separately. 



nary pulsars. These could be explained by the fact that so 
far we know only a few of such objects. 

On the other hand one can argue that before the sec- 
ond supernova explosion the stars are tidally locked, which 
would prevent the newly born neutron star from becoming 



a pulsar, provided that the orbit before the explosion is 
wide enough. There are two effects bracketing the size of 
orbit of binaries before second supernova explosion, and 
they have chance (depending on evolutionary parameters) 
to select binaries, which are tidaly locked before second su- 
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Fig. 3. Cumulative distribution of projected distances on the sky of compact object binary mergers. The left panels 

correspond to the case of propagation in the empty space, and the right panels to the propagation in the potential of 
a massive galaxy. The top panels show the case of no kick velocities cr„ = Okms~^, the middle panels correspond to 
(jy = 200kms~^, and the lower panels represent the case = 800kms~^. In each panel the solid line, long dashed , 
and short dashed lines represent the the distributions of NS-NS, BH-NS, and BH-BH binaries, respectively. 
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pernova explosion and are not too tight at that time. The 
first selection occurs during the first supernova explosion 
which tends to disrupt the widest binaries, as their bind- 
ing energy is small. After the first supernova explosion 
the population of binaries is reduced to these which are 
close, with relatively small orbits, and which are in general 
highly eccentric (the kick has random direction). As the 
system evolves toward second supernova, the secondary 
star enters consecutive giant stages and there is a good 
chance that it will fill a significant part of its Roche lobe 
and become tidaly synchronized. On the other hand the 
system can not be too tight to allow for the secondary gi- 
ant evolution and thus minimum separation is set by the 
size of evolved secondary. 

There could also be a gap in the distribution of com- 
pact object formed in supernova explosion, because of the 
fall-back of matter onto a compact object in the super- 
nova explosion (Fryer, 1999). If such a gap really exists 



than the number of BH-NS binaries is smaller than our 
simulations predict. 

The properties of the population of compact object 
binaries can be found from Figure ^ where we plot the 
distribution of different types of objects in the plane: v 
-the center of mass velocity, and tmrg - time to merge. 
In the the top panel, corresponding to the case of no kick 
velocities, = km s~^ , a correlation between v and tmrg 
is apparent. The center of mass velocities in this case are 



due to mass ejection from the system (Blaauw, 196C), and 
is of the same order of magnitude as the orbital velocity. 
The lifetime due to gravitational wave energy loss has been 



calculated by Peters (1964), and it scales like tmrg oc a 



where a is the orbital separation. The orbital separation is 
inversely proportional to the square of the orbital velocity, 
hence tmrg , which explains the trend in the top 

panel of Figure I With increasing the kick velocity only 
very tight and eccentric systems survive and their lifetime 
becomes shorter. The typical center of mass velocity is 
now determined by both the orbital velocity at the time 
of supernova explosion, and the kick velocity. This moves 
the center of the distribution towards shorter merger times 
and higher velocities with increasing the kick velocity. The 
distributions become wider and the correlation between v 
and tmrg vanishes for the case of high kick velocities. 



Similarly to our previous work ( [Bulik et al., 1999aD we 
present the distribution of compact object mergers pro- 
jected, distances on the sky for two cases: (i) propagation 
in the potential of a massive galaxy like the Milky Way, 
and (ii) propagation in empty space, corresponding to a 
small galaxy. The results the case (i) and (ii) are shown 
in the right and left panels of Figure |[ respectively. Gen- 
erally for the same width of the kick velocity distribution, 
the spatial distribution of more massive binaries is tighter 
around the origin for the case (ii). This is due to the fact 
that the merger time for the more massive objects is typi- 
cally smaller, see also Figure |[ The median projected dis- 
tance for the BH-BH binaries is about ten times smaller 



than that for the NS-NS binaries. Increasing the kick ve- 
locity leads to two effects: decrease of the lifetime of a 
compact object binary and increase of its center of mass 
velocity. The decrease of the lifetime, however, is more 
important in determining the spatial distribution of these 
objects. The median distance to the merger decreases sig- 
nificantly with increasing the kick velocity, however, a long 
tail of the distribution extending to large distances re- 
mains. 

In the case(i) (propagation in the potential of a mas- 
sive Galaxy) the increase of the width of the kick velocity 
distribution leads to division of the population into two 
groups: the bound and the escaping binaries. The spa- 
tial distribution of bound binaries follows that of the stel- 
lar matter in galaxy. This division is responsible for the 
breaks in the cumulative plots in the right panels of Fig- 
ure ^. There is a small difference between the populations 
of different types of compact object binaries, regardless 
of the kick velocity. With increasing the kick velocity the 
fraction of unbound binaries becomes larger. This fraction 
also depends on the type of binary, the escaping fraction 
of binaries containing black holes is smaller than that of 
the NS-NS binaries. For the case of the highest kick veloc- 
ity distribution width considered here, = 800kms~^ 
less than 10 percent of black hole neutron star binaries 
become unbound in comparison to almost 30 percent of 
the double neutron star systems. 

4. Discussion 



We have presented an extension of our previous work (Bu- 
[ik et al., 1999a) where we analyzed the distribution of 



neutron star binaries around galaxies. Here we have mod- 
ified the code to produce binaries which are heavier, e.g. 
contain black holes. The class defined as NS-NS binaries 



in this paper is broader than the binaries in ( Bulik et al. 



1999a| ), since it includes objects with masses up to SM©, 
while the heaviest compact objects in our previous paper 
had a mass of 2.4 Mq, and the bulk of them were 1.4 M© 
neutron stars. The trend, which we also noted before, that 
heavier binaries tend to stick closer to the host galaxies is 
evident when comparing the distribution around a mas- 
sive galaxy for the NS-NS class shown in this paper with 



the distributions shown in (Bulik ct al., 1999a). For the 
highest kick velocity distribution width considered here, 
(Jv = 800kms~^, the fraction of escaping NS-NS binaries 
is ~ 25% he re, in comparison to almost 40% in (Bulik 
let al., 1999a| 



When interpreting the results of Figure ^ one has to 
remember that for reasons stated above the kick veloc- 
ity probably decreases with increasing mass of the com- 
pact object. Thus, the distribution of NS-NS binaries 
lies somewhere between the case a^, = 200 km s"^ and 
(Jv ~ 800kms~^ (the middle and the lower panel in Fig- 
ure m and the distribution of BH-NS binaries must be 
somewhere between the cases = Okms~^ and cr„ = 
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200 km s-^ Wc find that BH-NS binaries (and also BH-BH 
systems) will merge within the host galaxies provided that 
they are massive. In the case of smaller galaxies, the BH- 
NS mergers will take place typically a few times closer to 
the host galaxy than the NS-NS mergers. However, there 
is still a long tail of the distribution, with a significant 
fraction extending to large distances. 

We have shown that the distances from the host galax- 
ies where compact object binaries merge decrease with in- 
creasing the mass of the binary. The reasons for such a be- 
havior are twofold. First, the lifetimes of heavier binaries 
are smaller, because of the increased gravitational wave 
energy loss. Second, the width of the kick velocity distri- 
bution in a supernova explosion is probably smaller when 
black holes are formed, which leads to smaller center of 
mass velocities. We conclude that BH-NS binary mergers 
are more likely progenitors of gamma-ray bursts than the 
NS-NS binaries. Their distribution closely follows that of 
the matter in massive galaxies, and also allows to explain 
the extreme energetics of some bursts. However, this con- 
clusion is so far based on a small sample of well observed 
GRB afterglows for long and hard bursts. 
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